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Abstract 

Mathematical diffraction theory is concerned with the analysis of the diffraction image 
of a given structure and the corresponding inverse problem of structure determination. 
In recent years, the understanding of systems with continuous and mixed spectra has 
improved considerably. Simultaneously, their relevance has grown in practice as well. In 
this context, the phenomenon of homometry shows various unexpected new facets. This 
is particularly so for systems with stochastic components. 

After the introduction to the mathematical tools, we briefly discuss pure point spectra, 
based on the Poisson summation formula for lattice Dirac combs. This provides an elegant 
approach to the diffraction formulas of infinite crystals and quasicrystals. We continue by 
considering classic deterministic examples with singular or absolutely continuous diffrac- 
tion spectra. In particular, we recall an isospectral family of structures with continuously 
varying entropy. We close with a summary of more recent results on the diffraction of 
dynamical systems of algebraic or stochastic origin. 
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1 Introduction 



Diffraction, originally observed and investigated in optics in the 17th century, developed into 
the main tool for structure determination of solids after X-rays were first used by von Laue and 
by Bragg in the early 20th century to determine the atomic structure of crystalline materials. 
Nowadays, a hundred years later. X-ray, electron and neutron diffraction continue to be among 
the most important techniques used in structure analysis (see [34J and references therein), 
complemented by direct imaging techniques such as electron and atomic force microscopy. 

The diffraction of a beam of X-rays, electrons or neutrons from a macroscopic piece of 
solid is a complicated physical process. The inverse problem of determining the structure from 
the diffraction intensities is even more involved, and rarely determines the structure uniquely. 
The presence of inelastic and multiple scattering, which is particularly prevalent in electron 
diffraction, makes it essentially impossible to arrive at a unified mathematical description of 
the process. In contrast, within the realm of kinematic diffraction in the Fraunhofer limit, 
there exist powerful mathematical tools to tackle the direct problem and to gain insight into 
the complexity of the corresponding inverse problem. 

While the theoretical description had been well developed for the case of ordinary (peri- 
odic) crystals and incommensurate phases, the discovery of quasicrystals in the early 1980s 
\86\ 154} l60l [68] not only required new mathematical methods, but also re-opened general 
questions concerning the possible manifestations of order and disorder in solids. In addition, 
improved experimental techniques make the diffuse part of the diffraction increasingly acces- 
sible, and a better understanding of the continuous diffraction is desirable, in particular in 
view of the implications on disorder. 

In this review, which is an extended version of lectures given by the authors at RIMS (Ky- 
oto) in summer 2010 [18], we summarise the mathematical development that was stimulated 
by the discovery of quasicrystals and has led to a more systematic approach to systems with 
all types of spectral components. 

2 Mathematical diffraction 

Since the pioneering paper by Hof |5T] , it has become evident that the formulation of math- 
ematical diffraction theory via measures is both powerful and versatile. As this is not the 
usual language of crystallography, we begin with a motivation and then summarise the math- 
ematical background for this review. 

2.1 Why measures? 

A mathematically satisfying formulation has to talk about infinite systems. Traditionally, they 
are often described by functions (representing the density of the scattering medium) or by 
tilings (whose decorations are meant to mimic the atomic positions). These two points of view 
are somewhat contradictory in the sense that each idealisation comes with its own (metric) 
topology (by which we mean a notion of 'closeness' of structures) which is not compatible 
with the other. 

To expand on this, almost periodic functions are compared via the supremum norm, and 
one function is close to another when each lives within an e-tube around the other, for some 
small e > 0. In particular, they are locally e-close everywhere. This distance concept makes 
no sense for the comparison of two Penrose tilings, say. There, one employs the local topology 
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instead, which says that two tihngs are close when, possibly after a translation of length at 
most e, they agree exactly on a ball of radius 1/e (while no condition is imposed in the 
complement of the ball). This, in turn, makes no sense for almost periodic functions, because 
they almost never agree exactly on finite regions. 

This dichotomy has led to rather fierce scientific disputes between the tiling school and 
the density function school, in particular in the years following the discovery of quasicrystals. 
Mathematically, however, this is a fight about nothing, because the two viewpoints can be 
reconciled by embedding them into a larger, and slightly more general frame. One possibility 
to do so is to introduce measures together with the so-called vague topology. Such measures 
comprise both almost periodic functions and tilings as special cases. Moreover, when re- 
stricting to one of these cases, the vague topology agrees with the other topologies mentioned 
before. 

We take this as sufficient motivation and justification to use a formulation with measures 
from the very beginning. In addition, measures provide the most natural frame for systems 
with disorder, as is well-known from probability theory. Meanwhile, the relevance of different 
topologies has been recognised in a more general setting; compare \27\ [7H [23] for more. Let 
us finally mention that a measure is a mathematical notion that is well suited to quantify 
both the distribution of matter and the distribution of (scattered) intensity in space, so that 
its appearance in our present (physically motivated) context is very natural. 

2.2 Mathematical preliminaries 

For simplicity, we introduce measures as linear functionals on continuous functions, and then 
connect them to the standard approach of regular Borel measures via the Riesz-Markov 
representation theorem; see |79j and references therein for background material. 

Let K. = Cc(M'^) be the space of complex- valued, continuous functions of compact support. 
A (complex) measure fx on is a linear functional on K, with values in M (or in C), with the 
extra condition that, for every compact set K C M'^, there is a constant such that 

1^(5)1 < Iblloo 

for all test functions g £ JC with support in K. Here, Halloo = ^^Px£k\9{^)\ is the supremum 
norm of g. If /i is a measure, the conjugate of fi is defined by the mapping g 1— )• fi{g). It is 
again a measure, denoted by fl. A measure is called real (or signed), when ft = /j,, and it is 
called positive when fj,{g) > for all g >0- For every measure //, there is a smallest positive 
measure, denoted by \fi\, such that \^J'{g)\ < IfJ-Hg) for all non-negative g £ JC. This is called 
the total variation (or absolute value) of fi. 

A measure ^ is bounded (or finite), if |/i|(l) = I^KM'^) (see below for the meaning of 
this notation) is finite. Otherwise, it is called unbounded, an example of which is given by 
Lebesgue measure A. The vector space of measures on M*^ is given the vague topology, which 
means that a sequence of measures {nn)nen converges vaguely to /i if lim„_>oo ^n(5) = l^-ig) 
in C for all g £ IC; compare [791 p. 114] for more on this topology in the context of linear 
functionals. Note that the term 'vague' has a very precise meaning this way, and vague 
convergence is a natural and powerful concept of analysis. The measures defined this way 
are in one-to-one correspondence with the regular Borel measures on M'^, by means of the 
Riesz-Markov representation theorem |79] . The cr-algebra of measurable sets is formed by the 
Borel sets, which is the smallest u-algebra that contains all open (and hence also all closed) 
subsets of M*^ in its standard topology. In view of this identification, we write ^{A) (measure 
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of a set A C M'^) and fj,{g) (measure of a function), and use the meaning interchangeably (for 
instance for characteristic functions of a set, where /^(l^) = fJ,{A)). 

The distribution of matter in Euchdean d-space is described by a measure u on M"^, where 
we assume an infinite system that is homogeneous and in equihbrium. In most cases, uj will 
be translation bounded, which means that, for any compact set K C W^, we have 

sup \ijj\{t + K) < oo . 



Moreover, we assume an amenability property of lo, namely the existence of its autocorrelation 
measure 

_ U}\„*UJ\„ 

-f = -f^ = UJ ® UJ := hm f , 

where Bji denotes the open ball of radius R around G M*^ and uj\j^ the restriction of uj to 
Br. Given a measure fi, its 'flipped-over' version Jl is defined via Jl{g) = nQj) for g £ IC, 
where g{x) = g{—x). The volume-averaged (or Eberlein) convolution ® is needed because u 
itself is an unbounded measure, so the direct convolution is not defined. For instance, if A 
denotes the standard Lebesgue measure (for volume), A * A is not defined, while A ® A = A. 
Note that different measures uj can share the same autocorrelation 7. This phenomenon is 
called homometry, and we shall see explicit examples later on. 

By construction, the measure 7 is positive definite (or of positive type), which means 
7(5 * 5 ) ^ for all g £ IC. As a consequence, 7 is Fourier transformable by general results 
|31j . The Fourier transform 7 exists and is a positive measure, called the diffraction measure 
of UJ. It describes the outcome of kinematic diffraction by uj in the sense that 7 quantifies how 
much scattering intensity reaches a given volume in d-space. By the Lebesgue decomposition 
theorem, relative to Lebesgue measure A, there is a unique splitting 

7 = 7pp+7sc + 7ac 

of 7 into its pure point part 7pp (the Bragg peaks, of which there are at most countably 
many) , its absolutely continuous part 73^ (the diffuse scattering with locally integrable density 
relative to A) and its singular continuous part 75^ (which is whatever remains). The last 
contribution, if present, is described by a measure that gives no weight to single points, but 
is still concentrated to an (uncountable!) set of zero Lebesgue measure. For a proof and 
background material, we refer to [791 Sec. 1.4]; see also [7] and references therein. 

Systems with 7 = are called pure point diffractive. Important examples are perfect 
crystals and quasicrystals [63], such as the icosahedrally symmetric AIMnPd alloy that pro- 
duces the diffraction image of Figure [TJ Mathematical examples of all spectral types will be 
discussed below. For general background, we refer to [26], and to [12] in particular. The 
increasing need for a better understanding of diffuse scattering is evident from [9T] and ref- 
erences therein. 

It is clear that the diffraction measure is a unique attribute of the structure described by 
UJ (under the mild assumption that 7 exists, which is a realistic assumption from the physical 
point of view, and very natural mathematically). In contrast, the inverse problem of deter- 
mining UJ from 7 is generally non-unique and a hard problem to solve, both mathematically 
and practically. 
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Figure 1: Experimental diffraction pattern of an icosahedral AIMnPd quasicrystal along the 
fivefold direction (courtesy C. Beeli). 



2.3 Lattice Dirac combs 

For simplicity, we will, whenever possible, explain the different scenarios with (weighted) 
Dirac combs on Z. This means that, given a bi-infinite sequence w = (u^n)„gZ' consider 



where 5n is the normalised Dirac measure located at n. (Recall that 5x, for a fixed a; G M, 
is defined by Sx{g) = g{x) for arbitrary g e K.). The weights w{n) = Wn are assumed to be 
bounded, and we use the notations w{n) and Wn interchangeably. A simple calculation shows 
that the corresponding autocorrelation (which will indeed exist in all examples discussed here) 
is then of the form 




(1) 



neZ 




(2) 
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with the autocorrelation coefficients 

I N 

r/(m) = lim V w{n)w{n-m). (3) 

n=-Af 

In hne with our previous remark, we wih assume that the hmit exists for all m G Z, which 
is equivalent to the existence of 7 in this case. A stochastic setting will actually require the 
existence in a suitable probabilistic sense; we postpone this point of view until Section [5l 

Let us now illustrate the possible spectral types by means of important and characteristic 
examples. 



3 Pure point spectra 

The Dirac comb of a general point set S is defined as 5g := Ylixi^s with the normalised 
point (or Dirac) measure at x. If 7^ C M is a lattice (meaning a discrete subgroup of such 
that the factor group / F is compact), the corresponding Dirac comb 5r itself is Fourier 
transformable via the Poisson summation formula (PSF) 

6p = dens(r) 5^. , (4) 

where F* = {x £ M.'^ \ xy ^ F for all y S F} denotes the dual lattice of F and dens(-r) the 
density of F; see [33\ [27] and references therein for details. This formula is fundamental in 
many ways (and disciplines), and can easily be proved in two steps; see f7] for an explicit 
exposition. Let us explain what ([!]) means for crystals and quasicrystals. 



3.1 Crystallographic systems 

A perfect (infinite) crystal with F as its lattice of periods can be described by the measure 
oj = fi * 5p, where ^ is a suitable finite measure. A convenient (though not unique) choice 
for /i is the restriction of a; to a fundamental domain of 7^. A simple calculation leads to the 
autocorrelation 

7 = dens(-r) {fi * Jl) * 5p 
because 5p = 5p and 5p ® 6p = dens{F) 5p. The Fourier transform of 7 exists and reads 

7= (dens(r'))^ |/Up (5^. 

by an application of the convolution theorem together with the PSF (jl]). Note that /I is a 
uniformly continuous and bounded function that is evaluated only at points of the dual lattice 
F* . Different admissible choices for the measure fi lead to different such functions that agree 
on all points of F* , so that the result does not depend on the choice made. The measure 7 is 
a pure point measure, as one expects for lattice periodic measures uj. 



3.2 Model sets 

Lattice periodic point sets form a special case of the larger class of regular model sets [731183]. 
which also lead to pure point diffraction measures. One of the simplest non-periodic examples 
in one dimension emerges from the silver mean substitution rule 

a I— )• aba 
bh^ a 
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Figure 2: The silver mean point set A of Equation ([5]) as a projection of a strip-shaped subset 
of the lattice F = {{x,x') \ x ^ Z[-v/2 ]}. The shading marks the left endpoints of a-type (light 
grey) and 6-type (dark grey) intervals. 



which has inflation multiplier cr = 1 + \/2 . The latter is known as the silver mean and is 
the Perron- Frobenius eigenvalue of the corresponding substitution matrix (f g)- Note that a 
is a Pisot-Vijayaraghavan (PV) number, which means that it is an algebraic integer a > 1 
whose algebraic conjugates (apart from a itself) all lie inside the open unit disk; see [32] for 
background material and details. 

The natural geometric realisation of this system, starting from a bi-infinite fixed point of 
Q with legal seed a\a, is built via two intervals of length ratio a. If a represents an interval of 
length a and b one of length 1, their left endpoints constitute the silver mean point set 

A={xGZ[V2]\x' (5) 

the proof of which is not entirely trivial; it is spelled out in detail in Here, ' denotes 

algebraic conjugation in the quadratic field Q{y/2), as defined by y/2 i->- — \/2 . Equation ([5]) 
has a distinctive geometric meaning which is illustrated in Figure [21 

The algebraic aspects are encoded in a cut and project scheme (CPS), which we state for 
the general Euclidean follows: 



ad , Trad 



X 



U U U de 

1-1 



n{C) ^ C > (6) 

L y L* 



Here, £ is a lattice in M +™ with certain properties that are expressed via the images under 
the canonical projections vr and ttj^j. In particular, L = it{C) is a bijective image of £, while 
L* = 7rjjj^.(/^) is dense in internal space W^. Due to these properties, the *-map x i— )• x* 
is well-defined on L; see [73j for more. In the silver mean example, we have a CPS with 
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d = m = l, L = L* = Z[V2], and the T*r-map is given by algebraic conjugation, as mentioned 
above. 

In general, a model set for a given CPS is a set of the form 

A = {xe L\x* GW} (7) 

where is a relatively compact subset of M™; see Equation ([5]) for the silver mean case. A 
model set A is regular when the boundary dW of the window W has zero Lebesgue measure. 
The entire setting generalises, without significant complications, to locally compact Abelian 
groups as internal spaces [72} 1731 [83]. We will refer to this freedom later on, where the internal 
space will be based on the 2-adic numbers. 

The Dirac combs of regular model sets are pure point diffractive |51l [27] . This is a 
substantial theorem for which three different types of proofs are known. The most common one 
is based on the connection to dynamical systems theory |51 [ I83 [ [67]. another on a reformulation 
via almost periodic measures; see [STliZSKHTj and [ISl Sec. 4], and [TSl Lemma 6.6] for a related 
result in the context of subshifts. An even simpler approach follows a suggestion by Lagarias 
and is based on the PSF for the embedding lattice together with Weyl's lemma on uniform 
distribution [19]. The theorem is also constructive in the sense that it provides an explicit 
and computable formula for the diffraction measure of (5^, namely 

with Fourier module L® = 7r(£*) and amplitudes 



where lyy is the characteristic function of the window W. This formula has several gen- 
eralisations [831 [271 il9] - in particular to certain weighted Dirac combs, which we omit for 
simplicity. 

Model sets are widely used to describe and analyse diffraction images such as that shown 
in Figure [TJ Although real world quasicrystals will usually not be pure point diffractive, 
their average structure is well captured by this approach. In this regard, quasicrystals behave 
pretty much like ordinary crystals. 

Let us expand on the diffraction formula for the vertex set of the planar Ammann-Beenker 
(or octagonal) tiling [3]. This point set is a regular model set with L = where = 

exp(27ri/8) is a primitive eighth root of unity. The *-map (illustrated in Figure [3]) is defined 
by ^ I— >• which is an automorphism of the cyclotomic field Q{^,), so that L* = L. The 
lattice £ = {{x,x*) \ x £ is the Minkowski embedding of L, which is a scaled copy of 

in this case. The standard window is a regular octagon O of edge length 1, centred at the 
origin. The model set construction produces the point set 

Aab = {x G Zl + + + \x* eO} 

and its ★ -image A\-q of Figure HI The corresponding tiling emerges by connecting all vertices 
of distance 1 in ^ab- It has two prototiles, a square and a rhombus. 

The diffraction measure is calculated via the Fourier transform of the characteristic func- 
tion 1q . This leads to a dense (but countable) set of Bragg peaks whose intensities are locally 
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Figure 3: Action of the Z-linear ★-map on the generating elements of the eightfold module 





Figure 4: Patch of the Ammann-Beenker tiling and point set yl^g (left) and its lift to internal 
space (right). 
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Figure 5: Computed diffraction image of the Ammann-Beenker point set ^^B' ^^^^ 
details. 

summable. This means that, in any compact region, there are only finitely many peaks above 
any given positive threshold. A precise explicit calculation can easily be performed by one of 
the standard computer algebra packages. Figure [5] shows the result for a central patch of the 
Ammann-Beenker diffraction, with cutoff at 1/1000 of the central intensity. Here, a Bragg 
peak is represented by a small disk whose area is proportional to the intensity and whose 
centre is the position of the Bragg peak. The locations of the peaks above the threshold is a 
Delone set (in fact, it is even a Meyer set); see |88] for a general result in this direction. 

A related and very interesting question concerns the spectral type of systems that are 
defined by substitutions; see |78] for background. More recently, practically useful tests for 
pure pointedness of substitution systems have been developed, compare [2] and references 
given there, but the famous Pisot substitution conjecture (which states that the dynamical 
system associated to an irreducible Pisot substitution has pure point spectrum) and its higher- 
dimensional generalisations still remain a mystery. We will not discuss this point of view in 
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what follows, even though many of our examples will be defined by substitution rules. 

Let us close this paragraph with a brief general comment. A translation bounded measure 
u! also defines a dynamical system under the translation action of W^; see \H'6\ [22} [25] for 
background. If to is pure point diffractive, the corresponding dynamical spectrum is pure 
point as well (the converse also being true). This equivalence is well understood by now 
[IQIIMIIMIIIIIIIHIET] (a related result was already obtained in [THl Prop. 4.14]), but does not 
extend to general systems with continuous spectral components |42[ l9]. as we shall see later 
on. 

3.3 Homometry 

As mentioned previously, different measures may possess the same autocorrelation and hence 
the same diffraction. This phenomenon is called homometry ^7j. Clearly, St*u} (with t E M*^) 
as well as uj have the same diffraction as w, but non-uniqueness is generally not exhausted by 
this. Let us illustrate how it appears already among pure point diffractive systems. 

The simplest situation emerges for periodic Dirac combs on Z with rational weights. As 
an example, Griinbaum and Moore jSOl constructed homometric Dirac combs of the form 

5 

^ = Sgz*Y^ Cj Sj (8) 

j=0 

with positive integer weights Cj, which are thus 6-periodic. The two choices of Table [1] lead 
to the same autocorrelation. Even worse, these two cases have the same correlation functions 
up to 5th order, and differ only in higher orders. Nevertheless, the two Dirac combs are 
substantially different. Note that the diffraction measure, which is supported on Z/6, shows 
systematic extinction (the intensity vanishes on all points of the form k = 1/6 with £ = 2, 
3 or 4 mod 6). Such extinctions are an indication for non-trivial homometry of pure point 
diffractive systems. 

Table 1: Integer weights cj for the two homometric Dirac combs built from Equation ([8]); see 
[50] for details. 
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Let us go one step beyond and show a pair of homometric model sets. We use the CPS of 
the Ammann-Beenker tiling from above, but replace the octagonal window by one of the two 
polyominoes shown in Figure[6l where we follow [13] . The Dirac comb (5^ of the corresponding 
model set A with window W has an autocorrelation of the form 

1A= Yl ^(^)'^-' (9) 
zeA-A 

where yl — 71 is locally finite because yl is a model set. The autocorrelation coefficients in ([9]) 
are given by 

voUwniw - z*)) 

r]{z) = dens(yl) -)— ^ = dens{C) cvg^iz*) . (10) 
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Figure 6: Two polyominoes (left) and their common covariogram (right). The small dots 
mark the origin, see text for further details. 

Here, cvg,^ is the covariogram of W, defined by 

cvgvK(2;) = vol{Wn{x + W)) = {lw*l_w){x). 

This function is symmetric under refiection in the origin, and one also has the relations 
cvgj_,_y(/ = cvg^y for arbitrary t G M'^ as well as cvg^_p^^ = cvg^y. The covariogram for our 
two polyominoes is illustrated as a contour plot in Figure [H The polyominoes in the same 
figure are shown with a shifted overlay structure, which can be used by the reader to check 
the claimed homometry (for the displayed shift) on the basis of Equations ([9]) and (|10p . 

The two model sets constructed this way differ on positions of positive density. Depending 
on the length scale of the windows, they may or may not be locally equivalent via a mutual 
local derivation (MLD) rule |5j, but they are always homometric. Further details are discussed 
in \13\ I49j. If one has access to correlation functions of higher order, a distinction is possible. 
In our example, the 3-point correlations tell the two examples apart. Some rather general 
results in this direction were recently derived in \39\ l66] . 

4 Singular continuous spectra 

The probably best known singular continuous measure is the one that emerges from the 
middle-thirds Cantor set construction. Its distribution function F is shown in Figure [71 which 
is widely known as the Devil's staircase. This function is continuous and non-decreasing, but 
constant almost everywhere. More precisely, the underlying measure fi = dF is concentrated 
on the Cantor set C, which is an uncountable set of zero Lebesgue measure. The (positive) 
measure fi is singular continuous, with /i({x}) = for all x £ [0, 1] and ^(C) = 1. 

4.1 Thue-Morse sequence 

Let us now discuss a classic example from the theory of substitution systems that leads to a 
singular continuous diffraction measure with rather different features in comparison with the 
Cantor measure. This example has a long history, which can be extracted from |90 1 169 1 155 1 IT]. 
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Figure 7: Illustration of the construction of the middle-thirds Cantor set C (inlay) and the 
distribution function F for the corresponding probability measure on C. 



We confine ourselves to a brief summary of the results, and refer to [14] and references therein 
for proofs and details. 

The classic Thue-Morse (TM) sequence can be defined via the one-sided fixed point v = 
VqViV2 ■ ■ ■ (with Vq = 1) of the primitive substitution rule 

1 11 
i^ll 

on the binary alphabet {1, 1}. The fixed point is the limit (in the obvious product topology) 
of the iteration sequence 

1^11^ 1111 h-^ liilllli ^ ... — yv = g{v) = VQV1V2V3 . . . 

and has a number of distinctive properties [U [78] , for instance 

^ y, — ^ j^^sum of the binary digits of i 

• = and v^i+i = 

• V = VqV2V^ . . . and v = v^v^v^ . . . 

• V is (strongly) cube-free. 

Here, we define 1 = 1 and identify 1 with —1 for the later calculations with Dirac combs. A 
two-sided sequence w can be defined by 



Wi 



Vi, for i > 0, 
V-i-i, for i < 0, 
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which is a fixed point of , because the seed w_i\wq = 1|1 is a legal word (it occurs in 
and w = Q^iw). The (discrete) hull X = Xrpj^ of the TM substitution is the closure of the 
orbit of w under the shift action, which is compact. The orbit of any of its members is dense 
in X. We thus have a topological dynamical system (X, Z) that is minimal. When equipped 
with the standard Borel cj-algebra, the system admits a unique shift-invariant probability 
measure, so that the corresponding measure theoretic dynamical system is strictly ergodic 
[551178]. 

Any given u) G X is mapped to a signed Dirac comb (and hence to a translation bounded 
measure) w via 

w = Wn 5n ■ 

We inherit unique ergodicity, and thus obtain an autocorrelation of the form ([2]) with coeffi- 
cients rj{m) as in ([3]). Due to the nature of the fixed point an alternative way to express 
the coefficients is 



r]{m) = lim ^ Vn Vn+m 



n=0 

for m > together with r]{—rri) = r]{m). It is clear that r]{0) = 1, and the scaling relations of 
V lead to the recursions 

r]{2m) = r]{m) and 

1 (11) 
ri{2m+l) = --[r]{m) +r]{m + l)), 

which are valid for all m G Z. In particular, the second relation, used with m = 0, implies 
77(1) = — |, which can also be calculated directly. 

Since u is supported on Z, the diffraction measure 7 is 1-periodic, which follows from [Bl 
Thm. 1], but can also be seen from 7 = rjd^. Indeed, this implies h'j = 'j with the continuous 
function h{k) = e^'^^^ (which is 1 on integer arguments). Taking the Fourier transform and 
applying the convolution theorem backwards, one thus obtains 5i * 7 = 7. Moreover, 7 is of 
the form 7 = ^1*6^, where 

^ = 7|jQ together with r/(m) = J e^'^'"'^ d^(y) , 

the latter due to the Herglotz-Bochner theorem [561 Thm. 1.7.6]. One can now analyse the 
spectral type of 7 via that of the finite measure fi, where we follow [55]. Defining IJ{N) = 
Sm=-Af (^("^)) ' ^ two-step calculation with the recursion (fTT]l establishes the inequality 
E{AN) < |r(2Ar) for all N £ N. This implies Huin ^00 ^{N)/N = 0, wherefore Wiener's 
criterion [901 [56] tells us that /i is a continuous measure, so that 7 cannot have any pure 
point component. Note that the absense of the 'trivial' pure point component on Z is due to 
the use of balanced weights, in the sense that 1 and —1 occur equally frequent, and thus the 
average weight is zero. 

Let us now define the distribution function F by F{x) = fi[[0,x]) for x G [0, 1], which is 
a continuous function that defines a Riemann-Stieltjes measure [61:, Ch. X], so that dF = ji. 
The recursion relation for r/ now implies [55] the functional relations 



dF(f)±dF(^) = {_eosV)} 
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which have to be satisfied by the ac and sc parts of F separately, because /j,^^ _L /j,^^ in the 
measure theoretic sense; see [791 Thm. 1.20] or [611 Thm. VII. 2. 4]. Therefore, defining 

Jo 

we know that the coefficients r]^^{m) must satisfy the same recursions (|lip as i](m), possibly 
with a different initial condition rj^^^O). The Riemann-Lebesgue lemma [561 Thm. 1. 2. 8] states 
limm-s.±oo ??ac("^) = 0, with is only compatible with r/ac(0) = 0, because f/ac(l) = — |??ac(0) 
and ?7ac(2m) = ??ac(^) for all m € N, and hence ryac = 0. This means -Fac = by the Fourier 
uniqueness theorem, wherefore fx and hence 7 (neither of which is the zero measure) are 
purely singular continuous. The resulting distribution function is illustrated in Figure [H It 
was calculated by means of the uniformly converging Volterra iteration 

= 2 cos(7r?/))F„'(y) dy 

with Fq{x) = X. In contrast to the Devil's staircase, the TM function is strictly increasing, 
which means that there is no plateau (which would indicate a gap in the support of 7); see 
|14] and references therein for details and further properties of F. 

Despite the above result, the TM sequence is closely related to the period doubling se- 
quence, via the (continuous) block map 

(f): ll,llh^a, ll,llh^6, (12) 

which defines an exact 2-to-l surjection from the hull Xrpj^ to Xp^, where the latter is the 
hull of the period doubling substitution defined by 

Qp^ : a ab , h ^ aa . (13) 

Viewed as topological dynamical systems, this means that (Xpj,Z) is a factor of (Xrpj^,Z). 
Since both are strictly ergodic, this extends to the corresponding measure theoretic dynamical 
systems. The period doubling sequence can be described as a regular model set with a 2- 
adic internal space [281 EZ] and is thus pure point diffractive. This pairing also explains 
a phenomenon observed in [42j, as the missing part of the dynamical spectrum of the TM 
system is recovered via the diffraction measure of Xpd. 

4.2 Generalised Morse sequences 

The above example can be generalised to the family 

with k,^ G N, inspired by [57J. They define a class of systems which we will refer to as 
the gTM systems. All display purely singular continuous diffraction, which follows from 
completely analogous arguments [16]. The entire analysis is based on the structure of the 
autocorrelation, which reads 7 = r/ 5^ with r/(0) = 1 and the recursion relations 

r]{{k + l)m + r) = (a^^^ r/(m) + a^^^ ^^^.^ 7?(m + 1)) , 
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Figure 8: The strictly increasing distribution function F of the Thue-Morse diffraction mea- 
sure on the unit interval. 

with a^^j, = k + i — r — 2 min(A;, i,r,k + i — r). They are valid for all m G Z and < r < k + i. 
In particular, one has r](^{k + £)rnj = rj{m) for m G Z. 

Given G N, the distribution function F is defined by F{x) = 7([0,x]) for < x < 1, 
which extends to x G M via F{x + 1) = 1 + F[x). It is also skew-symmetric {F[—x) = —F{x)) 
and thus satisfies F(q) = q for all q G ^Z. The continuous function F possesses the uniformly 
converging series expansion [16j 

\ ^ rj(m) 

Fix) = x+ > smizvrmx) . 

^ mvr 

m>l 

Two further examples are shown in Figure [9l As mentioned above, F defines a purely singular 
continuous measure, for all values k,i €zN. 

Another analogy with the TM system is that the block map (j) of Equation (jl2p can still 
be used to induce a matching family of generalised period doubling sequences, which turn out 
to be conjugate to the Kronecker factor of the gTM system; compare [78j for the concept. 
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0.5 1 0.5 1 

Figure 9: The continuous and strictly increasing distribution functions of the generahsed 
Morse sequences with parameters (2, 1) (left) and (5, 1) (right). 

These sequences are defined by the primitive substitution rules 

that are all pure point diffractive. The latter claim can most easily be seen from the letter 
coincidence (in the sense of Dekking [37]) at the kth. position of the images. As before, this 
factor explores the pure point part of the dynamical spectrum of the gTM sequence. 

5 Absolutely continuous spectra 

The appearance of absolutely continuous diffraction spectra is usually seen as an indicator for 
randomness in the structure; see [23] for a discussion in the context of diffraction. Though 
this is perhaps generically true, there are also prominent deterministic sequences with such 
spectra, such as the Rudin-Shapiro sequence. 

In general, within the realm of random structures, one can only expect almost sure conver- 
gence results. In other words, most statements become measure theoretic in nature, though 
they are still completely rigorous. To tame randomness by means of strict laws is one of the 
tasks of probability theory, at which it is pretty successful; see [MIHS] for general background. 
Still, this point of view (which can be rather counterintuitive at times) is something one has 
to get used to. 

5.1 Coin tossing sequence 

The simplest example emerges from repeated coin tossing. Here, one obtains sequences w G 
{±1}^ (for instance with 1 for 'head' and —1 for 'tail') which may be considered as the 
outcome of an eternal coin tosser. In more modern (and slightly more general) terminology, 
one considers a family (^ra)^^^ of independent and identically distributed (i.i.d.) random 
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variables with values in {±1} and probabilities p (for 1) and 1—p (for —1). The ensemble 
of all possible realisations is Xg = {±1}^, which is equipped with a probability measure 
that emerges from the elementary probabilities via independence j44j . This gives a measure 
theoretic dynamical system that is called the Bernoulli shift [89]. It has (metric) entropy 
H{p) = -p\og{p) - (1-p) log(l-p). 

A random sequence W leads to a Dirac comb uj = WS^, which is now a translation 
bounded random measure with support Z. Its autocorrelation, if it exists, is of the form 
7b = Vb with 

VBim) ■■= lim WnWn+m ' = ' \ ]\ , JJ' 

N^<x, IJS+l y(2p-lY, m / 0. 

Here, the convergence (for m 7^ 0) is almost sure by the strong law of large numbers (SLLN) 
|43] . which means that this is the result for //^-almost all elements of {±1}^; see also [25] 
for further results on the basis of this type of argument. Note that the use of the SLLN can 
also be replaced by an application of Birkhoff 's ergodic theorem, because the Bernoulli shift 
is ergodic [SU]. The corresponding diffraction measure reads 

7^ ' = ' {2p-l)H^ + Ap{l-p)X, 

which follows from 7^ by an application of the PSF Q together with 5q = A. For the fair 
coin ip = this simplifies to 7g = A, which is thus our first example of a purely absolutely 
continuous diffraction measure. Note that the absence of a pure point component is a result of 
the (almost sure) balance between the weights 1 and —1, which is the probabilistic counterpart 
of the (deterministic) balance in the TM example. For a more general discussion of pure point 
spectra and entropy, we refer to [23] and references therein. 



5.2 Rudin-Shapiro sequence 

Let us contrast the coin tossing sequence with a deterministic example that derives from 
|821 185j . This sequence was originally constructed to show that the absence of pair corre- 
lations does not imply the presence of randomness. This has interesting consequences in 
diffraction theory, as pointed out in |53j . The modern formulation of the system is based on 
the substitution 

Q : ac, b ^ dc, ab, d ^ db . 

Since 5|a is a legal seed (it occurs in g'^{b)), one can construct a bi-infinite sequence u by the 
usual iteration procedure as 

b\a dbab\acab 1-^ . . . — > u = Q^{u) , 

where convergence is in the standard product topology. The hull (orbit closure) of u defines 
the quaternary Rudin-Shapiro system. Its reduction to a binary system is achieved by the 
mapping 

if : a,ci— s-l, 6, di— )-l, 

and the orbit closure of w := (p{u) = . . . 1111|1111 . . . defines the hull Xj^g of the binary Rudin- 
Shapiro system. The sequence w is illustrated in Figure [TOj For the equivalent description as 
a weighted Dirac comb on Z, we again use the identification of I with —1. 
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Figure 10: Illustration of the binary Rudin-Shapiro sequence. 



An alternative description of w uses the initial conditions w{—l) = —1, w{0) = 1 together 
with the recursion 

\{n), for £ e {0,1}, 

-l)"+^u;(n), for £ G {2,3}. 



The autocorrelation of the corresponding weighted Dirac comb Wj^g exists and turns out to 
be 7j^g = Sq. To prove this, one defines the coefficients 

?7(m) 1 __ w(n) w(n + m) I , \„ 

-dim)} JV^oo 2iV+l^f-'^ ^ ^ ^ ^\(-l)". 

An application of Birkhoff 's ergodic theorem to the quaternary Rudin-Shapiro system (which 
is strictly ergodic [78]) establishes the existence of all these limits. The recursion (jl4p now 
implies the recursion relations |16j 



r]{Am) = — 2 rjim) 



r/(4m + l) = ^ ri{m) + ^^-f^ ^{m) - \ 'd{m + l 



together with 



r/(4m + 2) = 0, 

r/(4m + 3) = 7?(m + l) - T?(m) + \§{m + l), 

^(4m) = 0, 

i?(4m+l) = ^~^4^^" ??(m) - ^^^j(m) + \d{m + l), 
t?(4m+2) = ■t:i^??(m) + i??(m+l), 

??(4m + 3) = r?(m + l) - §{m) + ^ t^(m + l). 

with initial conditions //(O) = 1 and 'i?(0) = 0. The set of equations closes in the sense that 
no new quantities occur. This was the reason for the introduction of the signed coefficients 
t?(m), which show up in the first block of equations. The system has the unique solution 
= together with r/(m) = 5m,o-, hence 7j^g = 5q. 
The diffraction measure for the binary Rudin-Shapiro system is thus given by 7j^g = A, 
which coincides with that of the coin tossing sequence for p = h; see also [53]. These two 



2' 

examples are thus homometric, despite the fact that one is deterministic (with entropy 0) 
while the other is fully stochastic (with entropy log(2)); see [151 [16] further details and 
discussions. 

5.3 Bernoullisation 

The homometry between the Dirac combs of the Rudin-Shapiro and the balanced coin tossing 
sequence raises the question how 'bad' the non-uniqueness of the inverse problem in this case 
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really is. A partial answer (to the negative) can be given by means of the 'Bernoullisation' 
procedure which was introduced in [i9l [T5] . 

Starting from a uniquely ergodic bi-infinite sequence S G {±1}^, its Dirac comb ujg = S 6^ 
possesses a unique autocorrelation 7^. Let us now consider the random Dirac comb 

W = ^ SnWn 5n , 

where {Wn)^^^ is once again an i.i.d. family of random variables with values in {±1} and 
probabilities p and 1—p. Another (slightly more complicated) application of the SLLN shows 
that Lo almost surely has the autocorrelation 

If 5* is the binary Rudin-Shapiro sequence, which is uniquely ergodic, a short calculation 
reveals that 

1 = ^0 = 7rs 

in this case, irrespective of the value of the parameter p £ [0, 1]. This way, we find a one- 
parameter family of homometric (or isospectral) structures whose entropy varies continuously 
between and log(2). The conclusion is that kinematic diffraction alone cannot distinguish 
order from disorder here [15] , 

6 Further directions 

After our brief discourse on the different spectral types of diffraction measures (by means of 
some paradigmatic examples), this section aims at indicating some more recent developments, 
which will again be explained informally by means of selected examples. 

6.1 Ledrappier's model 

Let us consider a prominent planar example of algebraic origin, which is due to Ledrappier 
[62] . It is defined as 

Xl = G {±1}^^ I W^W^+e-^Wx+e^ = 1 for all X G I?], 

where ei and 62 denote the standard Euclidean basis vectors in the plane. Xl is a closed 
subset of the full shift {±1}^ and hence compact. It is also an Abelian group (under point- 
wise multiplication in our formulation, which follows [29]). As a dynamical system, it is thus 
equipped with the corresponding Haar measure /i^i which is positive and normalised so that 
//l(Xl) = 1. Obviously, the system has no entropy, because the knowledge of a configuration 
along one horizontal line determines everything above it. However, it is clearly not deter- 
ministic. In fact, essentially along any given lattice direction, it looks like a one-dimensional 
Bernoulli system [62] . It is thus said to have rank 1 entropy, which essentially means that 
the number of circular patches of a given size grows exponentially with its diameter, but not 
with its area. 

Given an element w G Xl, the corresponding Dirac comb 
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possesses /iL-almost surely the autocorrelation 7 and the diffraction measure 7 given by [29] 

7 = (Jg and 7 = A . 

The system is thus homometric with the two-dimensional Bernoulli system with p = ^ (coin 
tossing on Z^), and also with the direct product of two binary Rudin-Shapiro sequences. The 
similarity with the Bernoulli system goes a lot further, in the sense that also other correlation 
functions agree, although the systems differ for certain 3-point correlations; see [29] for details 
and extensions. 

This system is meant to indicate that higher-dimensional symbolic dynamics is good for 
a surprise, as is well known from [HI]. It is thus clear that the inverse problem becomes 
more complicated with growing dimension, as new phenomena show up. Another famous 
example is the (x2, x3) dynamical system, which shares almost all correlation functions with 
a Bernoulli system with continuous (S^ ) degree of freedom [29] . 

6.2 Random dimers on the integers 

Back to one dimension, let us briefly describe a system that was recently suggested by van 
Enter [9]. First, consider Z as a close-packed arrangement of 'dimers' (pairs of neighbours), 
hence without gaps or overlaps. There are two possibilities to do so. Next, give each pair 
a random orientation by decorating it with either (+, — ) or (— ,-|-), with equal probability. 
Identifying it with ±1, this defines the closed (and hence compact) set 

X = {w£ {±1}^ I M{w) C 2Z or M{w) C 2Z + l} , 

where M{w) := {i €z Z \ Wi = tfj+i}. Note that M{w) can be empty, which happens precisely 
for the two periodic sequences . . . H — | H — ... and ... — 1-| — \- . . .. One has an invariant 
measure on X that emerges from the stochastic process via the probability of the possible 
finite patches (which define the generating cylinder sets as usual). 

Turning a configuration w ^X into a signed Dirac comb with weights Wi G {±1} as before, 
another exercise with the SLLN shows that its autocorrelation almost surely exists. It is not 
difficult to derive [9] that 7 = (5o — ^(^i + 5_i), so that the corresponding diffraction measure 
is 

= (1 — cos(27r/c)) A . 

This is another example of a purely absolutely continuous diffraction measure. The Radon- 
Nikodym density relative to A is written as a function of A;. In contrast, the dynamical 
spectrum of this system contains eigenvalues, wherefore this is an analogue of the Thue- Morse 
system |42J, this time in the presence of absolutely continuous spectra. 

The difference can be rectified by a block map that is very similar to the map (j) encountered 
in ()12p . Defining Ui = —WiWi-^-i for i G Z maps w to a new sequence u, which almost surely 
has the diffraction 

7n = ldz/2 + 2^ 

of mixed type [9]. In particular, it displays the entire point part of the original dynamical 
spectrum at once, as does the pure point diffraction of the period doubling system [28^ [T7] 
for the Thue-Morse case |16j . 
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6.3 Renewal processes 

A large and interesting class of processes in one dimension can be described as a renewal 
process [H |2T]. Here, one starts from a probability measure on (the positive real line) 
and considers a machine that moves at constant speed along the real line and drops a point 
on the line with a waiting time that is distributed according to /i. Whenever this happens, 
the internal clock is reset and the process resumes. Let us (for simplicity) assume that both 
the velocity of the machine and the expectation value of /i are 1, so that we end up with 
realisations that are, almost surely, point sets in M of density 1 (after we let the starting point 
of the machine move to — oo, say). 

Clearly, the resulting process is stationary and can thus be analysed by considering all 
realisations which contain the origin. Moreover, there is a clear (distributional) symmetry 
around the origin, so that we can determine the corresponding autocorrelation 7 of almost 
all realisations from studying what happens to the right of 0. Indeed, if we want to know the 
frequency per unit length of the occurrence of two points at distance x (or the corresponding 
density), we need to sum the contributions that x is the first point after 0, the second point, 
the third, and so on. In other words, we almost surely obtain the autocorrelation 

7 = 5q + v + V (15) 

with u = fi + fj.*fi + fi*n*fj, + ..., where the proper convergence of the sum of iterated 
convolutions follows from [U Lemma 4]. Note that the point measure at simply reflects 
that the almost sure density of the resulting point set is 1. Indeed, v is a translation bounded 
positive measure, and satisfies the renewal relations (see [451 Ch. XI.9] or [8, Prop. 1] for a 
proof) 

= n + fi*]^ and (1 — /I)z? = /2, (16) 

where is a uniformly continuous and bounded function on M. The second equation emerges 
from the first by Fourier transform, but has been rearranged to highlight the relevance of the 
set S = {k \ Ji{k) = 1} of singularities. In this setting, the measure 7 of (fTSj) is both positive 
and positive definite. 

Based on the structure of the support of the underlying probability measure /i, one can 
determine the diffraction of the renewal process as follows. Let /i be a probability measure 
on M-i- with mean 1, and assume that a moment of of order 1 + e exists for some e > (we 
refer to [8j for details on this condition). Then, the point sets obtained from the stationary 
renewal process based on jj, almost surely have a diffraction measure of the form 

7 = 7pp + (1 - h) A, 

where his a, locally integrable function on M that is continuous except for at most countably 
many points (namely those of the set S). It is given by 

_ 2{\m?-^<m)) 
^^"^ ~ \i-m\' 

see also [8] for further results on the basis of this type of argument. Moreover, the pure point 
part reads 

j 6q, if supp(//) is not a subset of a lattice, 
6j/L, otherwise. 



7, 



pp 
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where WL is the coarsest lattice that contains supp(/i). Proofs can be found in [8", 21]. 

The renewal process is a versatile method to produce interesting point sets on the line. 
These include random tilings with finitely many intervals (which are Delone sets) as well as 
the homogeneous Poisson process on the line (where [i is the exponential distribution with 
mean 1); see O Sec. 3] for explicit examples and applications. In particular, if one employs 
a suitably normalised version of the Gamma distribution, one can formulate a one-parameter 
family of renewal processes that continuously interpolates between the Poisson process (total 
positional randomness) and the lattice Z (perfect periodic order). 

6.4 Point sets from random matrices 

Another interesting class of random point sets derives from the (scaled) eigenvalue distribution 
of certain random matrix ensembles; see [21] and references therein. The global eigenvalue 
distribution of random orthogonal, unitary or symplectic matrix ensembles is known to asymp- 
totically follow the classic semi-circle law. More precisely, this law describes the eigenvalue 
distribution of the underlying ensembles of symmetric, Hermitian and symplectic matrices 
with Gaussian distributed entries. The corresponding random matrix ensembles are called 
GOE, GUE and GSE, with attached /3-parameters 1, 2 and 4, respectively. They permit an 
interpretation as a Coulomb gas, where /3 is the power in the central potential; see [U [7l] for 
general background and [41 ^ B6] for the results that are relevant here. 

For matrices of dimension N ^ the semi-circle has radius sJlN jix and area A^. Note that, in 
comparison with [71], we have rescaled the density by a factor Xj ^/l:^ so that we really have a 
semi-circle (and not a semi-ellipse). To study the local eigenvalue distribution for diffraction, 
we rescale the central region (between ±1, say) by \J2N jix. This leads, in the limit as — )• oo, 
to an ensemble of point sets on the line that can be interpreted as a stationary, ergodic point 
process of intensity 1; for /3 = 2, see [U Ch. 4.2] and references therein for details. Since 
the process is simple (meaning that, almost surely, no point is occupied twice), almost all 
realisations are point sets of density 1. 

It is possible to calculate the autocorrelation of these processes, on the basis of Dyson's 
correlation functions [41]. Though the latter originally apply to the circular ensembles, they 
have been adapted to the other ensembles by Mehta |71j. For all three ensembles mentioned 
above, this leads to an autocorrelation of the form 



where / is a locally integrable function that depends on /3; see [21] for the explicit formulas. 

The diffraction measure is the Fourier transform of 7, which has also been calculated in 
[m [7l]. Observing (5o = A and A = <5o, the result is always of the form 



where b = f. The Radon-Nikodym density h depends on /? and is summarised in [21j. 
Figure [TT] illustrates the result for the three ensembles. 

A similar approach is possible on the basis of the eigenvalues of general complex random 
matrices. This leads to the ensemble studied by Ginibre [71], which is also discussed in |21j . 
One common feature of the resulting point sets is the effectively repulsive behaviour of the 
points, which leads to the 'dip' around for 7. 




(17) 




(18) 
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Figure 11: Absolutely continuous part of the autocorrelation (left) and the diffraction (right) 
for the three random matrix derived point set ensembles on the line, with (3 £ {1,2,4}. On 
the left, the oscillatory behaviour increases with f3. On the right, f3 = 2 corresponds to the 
piecewise linear function with bends at and ±1, while /3 = 4 shows a locally integrable 
singularity at ±1. The latter reflects the slowly decaying oscillations on the left. 



6.5 Random clusters 

Let us continue by considering the influence of randomness on the diffraction of point sets 
and certain structures derived from them in Euclidean spaces of arbitrary dimension. Here, 
we start from a single point set A C W^, which is then randomly modified by replacing each 
point by a (possibly complex) finite random cluster. This situation is still manageable (via the 
SLLN) when A is sufficiently 'nice', for instance if A is of finite local complexity and possesses 
an autocorrelation, which is then of the form 7 = '^^^y^_y^r]{z) 6z- More generally, one can 
analyse this situation in the setting of stationary ergodic point processes [471 El EHl Ell [65] . 
which treats almost all realisations at once and permits a larger generality for the sets A, 
though the clusters will then be restricted to positive or signed measures. 

Given such a point set A, its (deterministic) Dirac comb 6y^ is turned into a random Dirac 
comb 

by means of the i.i.d. family of random measure {fix)x£A with common law Q and representing 
random variable f2. Here, we assume that the expectation EQ(|f2|) is a finite measure, and 
that Eq((|]7|(M'^))^) < c«. Under some mild (but somewhat technical) conditions |8|, one 
now obtains the autocorrelation 

^(O) (^^0 (]Eq(^) ^ i^)) * 7 + dens(^) (EQ(f2 * 17) - EQ{n) * E^)) * 60 

and hence the diffraction 

^(n) |Eq(1?)|2.^ + dens(yl)(EQ(|l?|2)-|EQ(/5)|2) -A. 

The diffraction of the modified structure emerges from the original one by a modulation of 
7 and the addition of an absolutely continuous contribution, which in essence is the Fourier 
transform of the covariance of the representing random cluster 17. 
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This approach comprises a wide range of models, including the random weight and the 
random displacement model [52] as well as decorations by random clusters. As mentioned 
above, the next level of generality replaces the deterministic set yl by a general ergodic point 
process |47] . This way, both the underlying set (core process) and the modification (cluster 
process) are described in terms of stochastic processes; see [SSJ [35] for general background 
and [8] and references therein for details. This point of view is a promising starting point for 
further investigations. 

7 Concluding remarks 

Our informal exposition was meant to demonstrate that mathematical diffraction theory 
provides useful tools for the analysis of deterministic and random systems, both for prac- 
tical applications in crystallography and materials science, and for theoretical questions in 
harmonic analysis and dynamical systems theory. While the majority of the crystallographic 
literature concentrates on the pure point case, we have shown that also continuous spectra are 
explicitly accessible in relevant cases, and certainly deserve more attention from this point of 
view. Merging methods from harmonic analysis and dynamical systems with well-established 
procedures from point process theory might be a good path to proceed. However, this does 
not only concern structures with some degree of disorder; there are also classes of completely 
deterministic systems whose diffraction is still not fully understood. A prominent example 
is provided by tilings with continuous symmetries such as the pinwheel tiling ^80j, which has 
circularly symmetric diffraction [811 [75] that resembles a powder diffraction pattern [TOl [XT]. 

The inverse problem of structure determination is already a formidable problem in the 
realm of pure point diffractive systems, due to the existence of non-trivially homometric 
structures. As we have shown above, the inverse problem becomes significantly more involved 
in the presence of disorder, including the potential insensitivity to quantities such as entropy. 
Although there is quite some knowledge about this problem in the point process community, it 
is fair to say that solutions to the inverse problem or satisfactory classifications of homometry 
classes are not in sight. 

The setting is by no means restricted to lattice systems, which were mainly chosen for 
ease of presentation and concreteness of results. Also extensions to higher dimensions are 
possible, where one has to expect new phenomena (such as the lower rank entropy) that 
further complicate the picture. Although the general theory of point processes is highly 
developed, the treatment of stochastic systems with interaction becomes difficult as soon as 
concrete results are desired. This is already the case for general equilibrium systems, as they 
require the full machinery of Gibbs measures [58[ [59| [70| [30[ |92j. Their analysis from the 
point of view of mathematical diffraction theory is still in its infancy, although examples such 
as the Ising lattice gas show (see [20] and references therein) that explicit results are possible. 
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